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Abstract 

A novel density functional, which accounts correctly for the equation of state, 
the static response function and the phonon-roton dispersion in bulk liquid 
helium, is used to predict static and dynamic properties of helium droplets. 
The static density profile is found to exhibit significant oscillations, which 
are accompanied by deviations of the evaporation energy from a liquid drop 
behaviour in the case of small droplets. The connection between such os- 
cillations and the structure of the static response function in the liquid is 
explicitly discussed. The energy and the wave function of excited states are 
then calculated in the framework of time dependent density functional the- 
ory. The new functional, which contains backflow-like effects, is expected to 
yield quantitatively correct predictions for the excitation spectrum also in the 
roton wave-length range. 
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I. INTRODUCTION 



The study of Helium droplets, as a prototype of finite quantum fiuid, has been the object 
of extensive experimental and theoretical investigations in the last decade (see Ref. for 
a recent review). The ground state and the excited states of pure ^He droplets, as well as 
of droplets doped with atomic and molecular impurities, have been studied using different 
theoretical approaches. However, the overlap between theories and experiments is not yet 
satisfactory. 

In this paper we investigate static and dynamic properties of ^He droplets using a density 
functional theory. We use a functional recently introduced by the Orsay-Trento collabora- 
tion p[. It is an extension of a previous phenomenological functional 0, which has been 
extensively used in the last years to study the static and the dynamic properties of inhomo- 
geneous phases of liquid helium (free surface, films, droplets). The Orsay-Trento functional 
improves, by construction, the description of relevant properties of the bulk liquid, such as 
the g-dependence of the static response function and the phonon-roton dispersion. We study 
here how this improvement affects the results for the static and dynamic properties of he- 
lium droplets. In particular, we look for deviations of the density profile and the evaporation 
energy from the liquid drop behaviour and we explore the effects of backflow correlations on 
the excitation spectrum at large momentum transfer. 

The paper is organized as follows. In Sect. || we present briefly the functional (for a 
complete discussion we refer to 0). In Sect. |T| we apply it to the ground state of the 
droplets. Then, in Sect. fV|, we sketch the formalism of time dependent density functional 
theory and present some results for excited states. 

II. DENSITY FUNCTIONAL 

A flrst systematic analysis of the ground state properties of ^He and ^He droplets in the 
framework of density functional theory was done in Ref. 0] , using a simple functional based 
on a zero-range effective interaction. The same functional was later generalized to include 
flnite-range effects 0. This functional has been applied to investigate static properties of 
several inhomogeneous systems, such as helium fllms and wetting phenomena |]5|, vortices 
in bulk liquid or droplets [0]. 

In density functional theory the energy of the system is written as 

E = jdvn[^,^*] , (1) 

with the complex function \i/ in the form 

v|/(r,t)=<l>(r,t)expQ%,t)) . (2) 

The real function $ is related to the diagonal one-body density by p = while the phase 
S flxes the velocity of the fluid through the relation v = {l/m)VS. The ground state of 
the system is found by minimizing the energy with respect to the density, while the time 
dependent theory is formally equivalent to the RPA equations with an effective particle-hole 
interaction, as we will see in Sect. We use the phenomenological density functional of 
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Ref . , containing a few parameters which are fixed to reproduce bulk properties of hquid 
hehum. It has the form: 



E = E^''^^ [p, v] + E(^) [p] + E^'f^ [p, v] , (3) 
where the first term is the kinetic energy of non interacting bosons, 

r h'^ 

E('=*")[p,v] = yc/r — |V^(r)|2 

= /rfr{|^-(VVp)^ + ^p(r)|v(r)r}, (4) 
the correlation energy is given by 

E(^)[p]= /rfr{i/rfr'p(r)V^(|r-r'|)p(r') + |p(r)(pr)^ + |p(r)(Pr)' 
- jdv' F{\v- r'\) [l - Vp(r) ■ Vp(rO 




(5) 

POs J \ POs J " 

and, finally, the backfiow energy IS 

E^'^\p,^r] = Jdrdr'Vj{\r-T'\) p{r)p{r') [v(r) - v(rO]' • (6) 

The first term in the kinetic energy, which depends on gradient of the density, is a quantum 
pressure; it corresponds to the zero temperature kinetic energy of non-interacting bosons of 
mass m. The two-body interaction V/ in the correlation energy E^'^^ is the Lennard- Jones 
interatomic potential, with the standard parameters a = 2.556 A and e = 10.22 K, screened 
at short distance (V^ = for r < h, with h = 2.1903A). The two terms with the parameters 
C2 = -2.411857 X 10^ K and C3 = 1.858496 x 10^ K account phenomenologically for 
short range correlations between atoms. The weighted density p is the average of p(r) over a 
sphere of radius h. Those terms are very similar to the functional of Ref. P]. The last term 
in E^'^\ depending on the gradient of the density in different points, is new and has been 
added in order to improve the description of the static response function in the roton region. 
The function F is a simple gaussian, F(r exp(— r V£2) with ^ = 1 A, while 

as = 54.31 A^ and pos = 0.04 A~^. The energy E^^-^"^ contains an effective current-current 
interaction accounting for backflow-like correlations. In Ref. the simple parametrization 

yj{r) = (711 + 712^^^) exp(-air^) + (721 + 722r^) exp(-a2'"^) , (7) 

was chosen in order to reproduce the phonon-roton dispersion in bulk liquid. The parameters 
are given in Table |. 

A detailed discussion about the meaning of the different contributions to functional @ 
and about the choice of the parameters is given in Ref. 0. Here we apply this functional 
to the statics and the dynamics of droplets. 
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III. GROUND STATE 



To find the ground state of the system one has to minimize the energy 



Eo = E(^\p] + jdv—iV^pf (8) 



with respect to the particle density p. This leads to the Hartree-type equation 



|-|^V^ + f/[p,r]|vpW = ■ (9) 

The Hartree potential is given by ?7[p, r] = SE^'^ySp{r). The chemical potential /i is in- 
troduced in order to ensure the proper normalization of the density to a fixed number of 
particles. 

For helium droplets the Hartree equation has to be written in spherical geometry. It 
becomes a unidimensional integro-differential equation which requires an iterative procedure 
to be solved. The resulting density profile, normalized to the bulk value po = 0.021836 A~^, 
for small droplets (10 < iV < 60) is shown in Fig. |I|. One clearly notice the appearance 
of oscillations both in the interior of the droplets and on the surface. The central density 
approaches the bulk density for large droplets, but it deviates significantly from the smooth 
behaviour predicted for a drop of compressible liquid @]. Such oscillations were not found 
with previous functional. 

Shell effects in ^He droplets, originating from Fermi statistics and shell closure in the 
single particle spectrum are expected to be important. On the contrary, the existence of a 
shell structure in bosons droplets is still under debate. Twenty years ago Regge and Rasetti 
suggested the possibihty to have "magic" droplets by looking at the structure of the 
elementary excitations in the bulk. In particular they showed that the surface of the liquid 
drop could act as a source of excitations producing static ripples on the density profile. In 
this approach the static ripples of the surface profile are connected to the pronounced peak of 
the static response function of bulk liquid at the roton wave length. The authors of Refs. 
also commented about the possibility that those oscillations could be washed out by the zero- 
point motion of the surface, whose thickness is expected to be larger than the interatomic 
distance. Indeed the subsequent theoretical predictions for the surface profile in helium 
droplets fljl0|-|12| , predicted a thickness of the order of 6 8 A and no density oscillations. 
A systematic analysis of small droplets with Green's Function Monte Carlo techniques was 



presented by Melzer and Zabolitzky |T3| in a work entitled "No magic numbers in neutral 
^He clusters" . However, the problem has been re-opened by the most recent Diffusion Monte 
Carlo calculations by Chin and Krotscheck |jl^ which show clear structures in the density 
profile of droplets with 20, 40, 70, and 112 atoms. 

In order to check the relation between the density oscillations in Fig. ^ and the behavior 
of the static response function x(g) we calculate the latter using our density functional 
theory. It can be evaluated by taking the second functional derivative of the energy through 
the relation 
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where V is the volume occupied by the system. A major advantage of hquid hehum is that 
x{q) is known experimentally since it is also equal to the inverse energy moment of the 
dynamic structure function S{q, u) which is measured in neutron scattering. In the upper 
part of Fig. ^ the experimental value of the static response function at zero pressure [15 



is 



compared with the predictions of three functionals: the one of the present work and the 
ones of Refs. [^3- The static response function is strongly g-dependent, showing a peak at 
the roton wave length. According to the theory of Regge and Rasetti this peak is also related 
to density oscillations in the profile of the free surface and of droplets. In the lower part of 
Fig. ^ we show the predictions for the density profile of the free surface with the three func- 
tionals. Indeed, while the three surface profile have a similar thickness, of the order of 6 to 8 
A, they have a different structure. Only the functional of the present work, which reproduce 
completely the peak of xil) the roton wave length, exhibits density oscillations. We note 
also that the same physical effect can be understood in terms of the static structure factor 
S{q), which has also a peak in the roton wavelength region. The relation between the quan- 
tities x(g) and S{q) in the framework of the present density functional theory is discussed in 
Ref. 0. It is worth mentioning that small density oscillations at the liquid- vapour interface 
have been recently predicted in classical fluids interacting through short ranged potentials 
TB| . In that case, the connection with the behavior of the radial distribution function, g{r), 



has been investigated in detail. For classical fluids the static response function and the 
radial distribution function are related through the Ornstein-Zernike equation. 

The oscillations of the free surface profile are very small; they are not expected to give 
rise to measurable effects. The surface energy is practically the same for the three profiles 
in Fig. (about 0.27 KA~^) and, moreover, the available experimental data on the surface 
reflectivity fl^ are not sensitive to such small oscillations in the profile. In the case of small 
droplets, however, the oscillating structure of the surface can combine with finite size effects, 
i.e., the tendency to form closed shell of atoms. This seems the case for the profile in Fig. |l|. 
We notice that the structure in the density of the droplets is not related to a "solid-like" 
behaviour; the oscillations originate at the surface, while a solid structures are not favoured 
in general by the presence of surfaces. 

The comparison between our results for the density profile and the ones of ab initio 
calculations is given in Fig. ^ for droplets with 20 and 70 particles. The solid line is the 
result of the present density functional calculation, while the DMC results of Ref. |n 



are 



represented by solid circles. The DMC data exhibit more pronounced oscillations, but 
possible spurious effects of metastable states, slowing down the convergence in the Monte 
Carlo algorithm, cannot be completely ruled out [|1^ . Recently the same authors have found 
oscillations in p(r) even with a variational approach based on the HNC approximation 



(dashed line). Even though the HNC method underestimates the central density, it predicts 
oscillations with amplitude and phase in remarkable agreement with the ones of density 
functional theory. An even better agreement is found in the most recent DMC calculations by 
Barnett and Whaley [jl9| (empty circles), where the statistical error is significantly reduced 
with respect to Ref. |T^ . 

Our predictions for the energy per particle are shown in Fig. ^ (solid line), together with 
the results of previous Monte Carlo and density functional calculations. The energy is a 
smooth function of the particle number N. Indeed ^He clusters behave "grosso modo" as 
liquid droplets and the energy can be easily fitted with the liquid drop formula 
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^ = a, + asN-^"" + a.N-^"" + aoN~^ , (11) 

where the volume coefficient a„ is the chemical potential in bulk, the surface coefficient is 
fixed by the free surface energy, while Oc and oq can be taken as fitting parameters. The 
energy calculated with the density functional differs from the liquid drop fit by less than 
0.02 K for all droplets with > 30. This seems to rule out apparently any shell effect in 
the energy systematics. However, one must keep in mind that the most relevant quantity 
determining the mass distribution of droplets is the evaporation energy [E{N — 1) — E{Ny\. 
The latter is not as smooth as the energy per particle. In Fig. ^ we show the evaporation 
energy predicted by the density functional (solid line) and the one obtained with the liquid 
drop formula (PD (dashed line). The difference between the two curves is also shown (dot- 
dashed line). When the difference is positive the droplets are more stable than it is predicted 
by the liquid drop formula. We note clear oscillations, having decreasing amplitude and 
increasing periodicity as a function of N . The same kind of oscillations appear in the 
central density of the droplets, as seen in Fig. |l|. Since the distance between two crests of 
the surface ripples is practically constant and the droplet radius goes approximately like 
A^^/^, the period of oscillations of the central density as a function of A^, as well as the one 
of the evaporation energy, increases as 

These results support quantitatively the original pictures proposed in Ref. 0, where 
the surface ripples have been discussed in terms of soft sphere closed packing. However 
the predicted deviations of the evaporation energy from the liquid drop behaviour are very 
small (less than 0.1 K). This makes their detection quite difficult since the temperature of 
the droplets in the experimental beams is presently expected to be about 0.4 K. 

In conclusion, we have shown some evidence for the occurrence of shell effects in neutral 
^He droplets. They are predicted by a density functional theory and physically originate by 
the strong peak in the static response function of liquid helium at the roton wave length. 
The shell effects are visible in the density profiles and the evaporation energy as a function 
of the number of particles in the droplet. Our results are in agreement with recent Monte 



Carlo simulations [0 and seem to confirm the previous predictions by Regge and Rasetti 



IV. EXCITED STATES 

A. Density-density response 

The first theoretical investigations of the elementary excitations of "^He droplets have been 
based on a generalized Feynman theory within the framework of microscopic Monte Carlo 
calculations of the ground state p!0|JT^j20|JI^ . A systematic discussion, using a Diffusion 
Monte Carlo algorithm, has been presented in [Q. In general, within this framework. 



it has been found that the dispersion with transferred momentum of the collective energies 
approaches the Feynman spectrum for bulk ^He as the cluster size increases. The convergence 
is very fast and already at small sizes (A^ ~ 70) there seems to be a well defined roton 
minimum. The Feynman spectrum is known however to overestimate by almost a factor 
two the experimental roton energy. 
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Since the Orsay-Trento functional reproduces the bulk phonon-roton spectrum, it is 
interesting to explore what it gives for the dynamics of helium droplets. In the following 
we study the behaviour of monopole (L = 0) and quadrupole (L = 2) excitations at low q, 
where the effect of the velocity dependent term turns out to be negligible, as well as the 
behaviour of monopole excitations at large q, where the spectrum is shown to approach the 
correct phonon-roton dispersion. We discuss the deviations from the results of a liquid drop 
model and the comparison with previous theories. 

The formalism of the time dependent density functional theory has been already intro- 
duced, for instance, in Refs. pT|J2^ J^ for applications to the dynamics of the free surface, 
of droplets, and films of helium on solid substrates. The same theory is here developed 
using the formalism of Green's Functions, allowing for a direct evaluation of the dynamic 
response function. We consider elementary excitations which are induced by an external 
field that couples to the particle density in the droplet. For sufficiently weak external fields 
the response can be treated linearly within the Random-Phase approximation. The response 
function to a transition operator Q{r) is defined as 

n{uj) = J dvdv' Q\v)G{y, r'; cj)Q(r) , (12) 

where G is the retarded Green's function, which is defined in terms of the ground state 
(|0 >), excited states {\n >) and their corresponding excitation energies Uno, and the creation 
operator ^/'^(r) (assuming time-reversal invariance of the matrix elements): 

G(y,y'-uj) = -Y.< 0|^^(r)^(r)|n >< n|^t(r')^(r')|0 > 



1 1 

+ 



Lu;„o -uj -ir] Uno + on + ir] 



(13) 



From Eq. (0) it is clear that the excitation energies are at the poles of G, and that they 
produce sharp peaks in the imaginary part of TZ{uj). 

Within the RPA the Green's function is calculated from the equation 

G'(r,r';^) = G'(°)(r,r';a;) + j dv^dv^ G^^^r, ri; a;)i;,(ri, r2)G'(r2, r'; ^) , (14) 

where G^^^ is the Green's function for a helium droplet with atoms within the single- 
particle model, i.e., expressed in terms of single particle wave functions (fi and energies 

Si-. 

G(°)(r,r» = -N^l{v)^,{v')Y.^,^{v)^l{v')[ ^ - + - \ ) , (15) 

and Vph is the residual (particle-hole) interaction. The sum in Eq. (^) extends, in principle, 
to all states ipn-, including those lying within the energy continuum (see Appendix). Once 
this is obtained, the RPA equation can be solved as a matrix equation in coordinate space 



and, finally, the response function can be obtained from Eq. (p!2D . 

In principle, the external field to be used in density-density response is the plane wave, 
with transferred momentum q, Q{r) = Y^iLi ^^^ ^ ■ However, with the multipole expansion 
indicated in the Appendix, one can calculate separately the response to each multipole Ql = 
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Sill iL(?^j)^Lo(^i)- In the limit of low q this external field reduces to Ql = J2iLi^t^Lo{ri)- 
Within the previous formulation of the response we will include exactly the coupling with 
the particle continuum. This may have some importance, especially for the response at high 
momentum transfer. 



B. Residual particle-hole Interaction 

The particle-hole interaction entering Eq. ([T^ ) is obtained, within density-functional 
theory, from the second variation of the energy functional. For functionals depending just 
on particle density p(r) and its gradients it may be obtained in the following way. 

Let the density be written, in general, in terms of the single particle basis fi{r) as 

Pi^) = Y.PiMir)Vjir) . (16) 

ij 

The second variation of the correlation energy Ec with respect to Pij's, taken at the ground 
state, provides the two-body interaction Vph- 

(5^ E'-"^ \p\ 

<tj\V,,ir„r,)\kl>=—-^. (17) 

OpikOpji 

The variation in this equation may be calculated as 



SpikSpji J \6p{r2)6p{ri) J ^ ^ 6pik 6pji 



and, consequently, we obtain the well known result 



VA^i,r2) = L ./'^^ J . (19) 



^6p{r2)5p{ri] 



g.s. 



The Fourier transform of this quantity enters the definition of the static response function 
x{ci) in bulk liquid, as in Eq. (|r^). From a practical viewpoint, the functional derivative 
implies the calculation of derivatives of the integrands in Eq. (^) with respect to p and Vp. 

The inclusion of the backfiow term in the functional requires a generalization of the 
previous scheme. This term introduces a dependence on the current density j, defined as 

= ^ E (V - V')p.,^r(r)^,(r')|_' • (20) 

Since the current density vanishes for the ground state, the only term contributing to the 
residual (backfiow) interaction will be 



g.s. 
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where the index v spans the three components, and a sum over v is assumed. From this 
expression one finds for the back-fiow part of the residual interaction 



•^'"'^'^'-i^dSr^l (V.-VO(V.-V,, (22) 




g.s. 



where the parenthesis is in fact independent of v and the gradients act only on the single- 
particle wave functions. This last expression shows explicitly how the current- dependent 
term yields a velocity- dependent residual interaction. Indeed, Eq. ( |22|) is quite similar to 
the Skyrme residual interaction in nuclear physics and we will use a similar technique 



to calculate the associated RPA response. Summarizing, we have two contributions to the 
residual interaction 



= V^^ + ) , (23) 



where V^j^^ is the residual interaction, obtained with Eq. (|T9|), associated to the functional 
without current terms, while V^j^^'' is the contribution from these current terms. 

C. Results 

In this section we present the results obtained for TZ{uj) in the formalism presented above. 
We use the technique developed in Ref. p3| , p4| and already applied to ^He droplets in the 



context of a contact (zero-range) effective interaction in Ref. [^|. The present formalism 



generalizes that of p2[ by including both finite range and back-fiow effects. We expect these 



to be important for the droplet response at high transferred momentum g, while for low q 



we expect to recover a surface-mode systematics similar to that of Ref. |]22[. The present 
functional reproduces by construction the dispersion of the elementary excitations in bulk 
liquid ^He [0, i.e., the phonon-roton curve up to g ~ 2.3 Consequently, it is particularly 
appropriate for the study of the response of droplets in this region of g's. 

1. Low q results 

Fig. ^ shows the imaginary part of TZ{u!) for the low q limit of the monopole and 
quadrupole external fields, which are Q = Y.i'^f and Q = J2i^i'^2o, respectively, for the 
droplet with N=112 atoms. The monopole corresponds to a breathing mode, while the 
quadrupole is a surface vibration [^. It is clearly seen how a very intense peak appears in 
the RPA result (continuous line), corresponding to a collective oscillation of the ^He atoms. 
The independent-particle result (dashed line), obtained by neglecting in Eq. (|l^ the particle- 
hole interaction, is in both cases more fragmented, revealing that the RPA correlations play 
a crucial role in the dynamics of the system. 

The effect of the current dependent part of the residual interaction, Vpl^'^\ on the 
monopole and quadrupole excitations in this range of q is completely negligible. We have 
explicitly checked that the monopole result in Fig. ^ is not affected at all by Vp^^\ This is 
consistent with the fact that in bulk liquid the lowest excitation for g — > is the phonon- 
mode, which is not affected by backflow-like correlations. As concern the L ^ modes of 
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the clusters, Chin and Krotscheck |T8| have recently shown that, at low g, they approach 
the dispersion of ripplons on the free surface. The latter, again, is found to be practically 
unaffected by backflow correlations 

The static induced densities Sp^^\r) for the same cluster are plotted in Fig. |^. This 
figure proves the behaviour mentioned before: the monopole induced density has a node and 
penetrates the interior of the cluster, while the quadrupole is almost completely localized 
at the droplet surface. An interesting result shown by this figure is the appearance of small 
oscillations in Sp^^\ in the inner part, both for the monopole and the quadrupole. This 
behaviour is connected to the oscillations displayed by the droplet equilibrium density p, 
and is due to the repulsive core of the effective interaction (see also Ref. [^). In fact, 
with the use of a contact interaction the small oscillations disappear both in p and Sp. 
Similar oscillations have also been found using other microscopic methods, like Monte Carlo 
calculations . 

The dependence of the collective energy with the number of atoms in the cluster is shown 
in Fig. |. The comparison with independent particle predictions shows that the role of RPA 
correlations increases with A^. The deviations of the RPA predictions from the Liquid 



Drop Model (LDM) results (long-dashed line, see p2[) are due to finite-size effects and are 



important up to A^ ~ 100 for surface vibrations and up to A^ ~ 500 for compression modes. 



The energies are in nice agreement with those of Ref. |^ , confirming the expectation that 
the response at g — > is not affected by the long range part of the atom-atom interaction. 
In the following section we show how this is completely different for the high q region of the 
response. 



2. High q results 

In this section we will focus on the L = response. In this case, the hole and particle 
states which contribute to the Green's function (|15D have spherical symmetry and this allows 
us to take only the radial part of the (V — V) in V^*'^"^'' , as a consequence, only the monopole 
term in the multipole expansion of the second functional derivative of Eq. (^) is needed. 
The situation is different for higher multipoles, since then the residual interaction ^^^"^^ 
involves different multipoles of the second functional derivative of Eq. (|2^). This greatly 
complicates the numerical calculation and we leave this case for a future investigation. In 
this work we concentrate on the monopole response to show the influence of the current 
term in the high-g response of droplets. 

Figure ^ shows the contour lines of the surface Im[7^(g, u)] for the A^ = 112 droplet. One 
clearly appreciates that the zone of higher response strength, where the peaks are located, 
resembles the dispersion curve of the bulk ^He elementary excitations. The existence of a 
roton minimum is clearly seen around q = 2 A~^. However, the energy of this excitation is 
quite affected by the current term Vp^'^K The upper panel is the result when V^j^'^^ is not 
included and provides a minimum energy at ~ 15 K. When this current term is included the 
minimum energy moves to ~ 9 K. Fig. ^ shows how this picture changes for different cluster 
sizes. For A^ = 20 the number of atoms is not enough to develop the roton minimum while 
for A^ = 728 the result follows quite closely the bulk one. Similar behaviour has been found 
using other microscopic methods []T2|,|T^ ; being based on Feynman-like wave functions, those 
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calculations yield a roton minimum at approximately twice the experimental roton gap in 
bulk liquid. 

V. CONCLUSIONS 

We have presented static and dynamic calculations for pure "^He droplets using a new 
density functional theory developed by the Orsay-Trento collaboration 0]. The density 
functional is written in such a way that relevant properties of the uniform liquid (equation 
of state, static response function, phonon-roton dispersion) are accurately reproduced. The 
theory is suitable to study properties of non uniform states of liquid helium. 

The static density profile and the energy of helium droplets have been calculated as 
a function of the number of atoms. We have found small density oscillations associated 
with regular deviations of the evaporation energy from the smooth liquid drop behaviour. 
We have discussed these oscillations in connection with the structure of the static response 
function, as suggested by a previous model of Regge and Rasetti Our results compare 
well with recent ab initio Monte Carlo calculations. 

We have also calculated the energy and transition densities of monopole and quadrupole 
excited states. The transition densities at g ^ display small oscillations similar to those 
of the ground state density. As in the case of the zero-range interaction, the solution of the 
RPA equations reveals the importance of the finite size effects and the long range correlations 
in the determination of the excitation energies of the compression and surface modes. In 
particular, the deviation from the LDM results are important up to ~ 500 for compression 
modes and up to A^ ~ 100 for surface vibrations. 

The analysis of the monopole mode at high momentum transfer shows the existence of 
a roton minimum around q 2 A~^ for clusters with a number of atoms A^ > 20. As in 
the uniform liquid, the backfiow correlations, which are included phenomenologically in the 
theory, yield a sizeable decrease in the energy of this minimum, from Ril5 K to ~9 K. 

ACKNOWLEDGMENTS 

This work was partially supported by European Community grant ERBCHRXCT920075; 
by INFN, gruppo collegato di Trento; and by DGICYT (Spain) grant PB92-0021-C02-02. 

APPENDIX A: RPA IN THE CONTINUUM 

The sum in Eq. ( [ISD can be written in a more compact way, using the particle Green's 
function Qp corresponding to the particle Hamiltonian Hp, in coordinate space: 

5(p(ri,r2;t^) =< ri|— -^^ |r2 > , (Al) 

rip — uj 

as 

G(°)(r,r';cu) = -N^l{v)^o{T')(g{Y,T'-eQ + uj + i7])+g{Y,v'-ea~uj~ir])) . (A2) 
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The particle Hamiltonian Hp is defined from Eq. (|]) 

Hp = -^V' + U\p,v], (A3) 

with p the ground state density. 

For spherical droplets each of the functions in Eq. can be expanded in multipoles as 

G(r,r') = Y.GL{ry)Yl,,{r)Y,M{r') , (A4) 

LM 



and, as a result, Eq. (14) also separates for different multipoles. In a symbolic notation 

Gz. = ^ + 4° . (A5) 

The particle Green's function for each partial wave is most easily calculated by means 
of the scattering solutions to the effective radial Hamiltonian Hp ^ at energy u: 

Hp^LUiir) = ujul ■ (A6) 

It is 



gL[r,r;uj) = . -, AT 

fi^ rr'W[UL,WL) 



where u and w are, respectively, the regular and irregular at the origin radial solutions of 
(|A6|), and W{ul,wl) is their Wronskian. In this way we can obtain G^^^ without having to 
impose any truncation of the continuum. 
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FIGURES 



FIG. 1. Density profile of ^He droplets for 10 < iV < 60, normalized to the bulk value. 

FIG. 2. Static response function of liquid helium (above) and free surface profile (below). 
Circles: experimental data jl^; lines: density functional calculations with functional of Ref. [Q] 
(dotted), Ref. d (dashed) and the one of Eq. (2) (sohd). 

FIG. 3. Density profile of two droplets. Solid line: present work; solid circles: DMC calcula- 
tions of Ref. [|l4l; empty circles: DMC calculations of Ref. [O]; dashed line: HNC calculations of 



Ref. |18|. 



FIG. 4. Energy per particle versus N. Solid line: present work; dotted line: Ref. Q; dashed 
line: results with functional of Ref. 0]; empty circles: Ref. |13]; solid circles: Ref. |14]. 



FIG. 5. Evaporation energy. Solid line: present work; dashed line: liquid drop formula. 
Dot-dashed line: deviation from the liquid drop formula (axis on the right). 

FIG. 6. Imaginary part of the response function at zero momentum transfer (g = 0) of the 
droplet with N = 112 atoms. Part (a) corresponds to the monopole case for which the external 
field is Q = J2i'''h (b) to the quadrupole, with Q = J2i^i'^2o{ri)- See text. 

FIG. 7. Induced densities for the droplet N = 112 in the static limit, oj = 0, corresponding to 
the results of Fig. |6|. We used arbitrary vertical scale. The dashed line shows the ground state 
density of the same droplet with the labeled vertical scale. See text. 

FIG. 8. Energy of the collective peak at g = as a function of size for the quadrupole and 



monopole. The short-dashed line shows, for comparison, the results of Ref. |22|. The long-dashed 
line corresponds to the LDM prediction. 



FIG. 9. Contour lines of the surface lm[TZ{q,uj)] for the monopole excitation of the droplet 

ribj 

ph 



N = 112. Part (a) is without including the current-dependent term V\ while (b) corresponds to 



the complete functional. See text. 

FIG. 10. Same as Fig. 9b, with the complete functional, for the droplets with = 20 (a) and 
N = 728 (b). 
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TABLES 



TABLE I. Values of the parameters used in Vj(r), see Eq. (0). 



711 


721 


712 


722 


ai 


02 


-19.7544 


-0.2395 


12.5616 A-2 


0.0312 A-2 


1.023 A-2 


0.14912 A-2 
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